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Abstract 

The strange and antistrange quark distributions of the nucleon are less constrained by exper- 
imental data than the non-strange quark sea. The combination of light quark sea distributions, 
A(a;) = d{x) + u{x) — s{x) — s{x), originates mainly from non-perturbative processes and can be 
calculated using non-perturbative models of the nucleon. We have calculated A(a;) using the me- 
son cloud model, which, when combined with the relatively well known non-strange light antiquark 
distributions obtained from global analysis of available experimental data, enables us to make new 
estimates of the total strange sea distributions of the nucleon and the strange sea suppression 
factor. 



PACS numbers: 14.20.Dh, 12.39.Ba 



I. INTRODUCTION 



The strange and antistrange quark distributions (hereafter referred to as strange sea 
distributions) of the nucleon are of great interest since these distributions are important for 
many processes in high-energy hadron coUisions. For example, a precise understanding of the 
cross-section for W production at the Large Hadron CoUider (LHC) depends on the strange 
sea distributions in the small x region. However, the strange sea distributions are not well 
determined compared with those for the non-strange, light quark sea. New measurements 
of charged kaon production in deep-inelastic scattering on the deuteron by the HERMES 
Collaboration pQ enabled them to present new determinations of the helicity averaged and 
helicity dependent parton distributions of the strange sea. Previous to this, the severest 
constraint on the strange sea distributions has come from di-muon production in neutrino 
and antineutrino deep inelastic scattering (DIS), mostly from the CCFR |2J and NuTeV [3] 
Collaborations. 

The standard method to extract the parton distribution functions (PDFs) of the nucleon 
is to assume a certain form for the various distributions at an initial scale Qq ~ 1 — 2 GeV^ 
and to perform a global fit to available experimental data for high energy scattering processes 
with the help of QCD evolution, which is known to next-to-next-to-leading order (NNLO) 
in the running QCD coupling constant. Many groups have been using this method to obtain 
the PDFs of the nucleon, see for example |1HTU]. An alternate approach to performing 
the global fit, mainly developed by the Dortmund group O HO], is to assume only valance- 
like distributions for the valance quarks, gluons and sea quarks at a somewhat smaller 
starting scale of = 0.2 ~ 0.5 GeV^. The PDFs at high are then mainly dynamically 
(radiatively) generated from these inputs via evolution. In all these global fits an initial 
difference between the strange quark sea and the light quark sea (i.e. SU(3) breaking among 
the PDFs of the nucleon) is assumed. In the standard approach a strangeness suppression 
factor 

s{x) + sx 

= J, . , . (1) 
a[x) + u[x) 

is introduced, and in the Dortmund group's approach the strange sea is assumed to vanish 
at the starting scale. While this initial difference between the strange and light sea improves 
the of thfi fits, we note that the uncertainties in the final strange and antistrange distri- 
butions are much larger than those of the light quark sea. There is also very little in the 



theoretical literature to explain this initial difference. While it is generally believed that the 
mass difference between the strange and light quarks is responsible for this difference, the 
mechanism for generating the difference has not been clearly elucidated. 

In the last few years, there has been an extensive effort to fully analyze the available 
data in order to learn more about the strange sea. The NuTeV data has been analyzed 
at next-to-leading order (NLO) [11] and the resulting strange and anti-strange distributions 
reported. A global fit to data from the neutrino(antineutrino)-nucleon DIS, inclusive charged 
lepon-nucleon DIS and Drell-Yan experiments, with an emphasis on extracting the strange 
sea distributions, was performed by Alekhin, Kulagin and Petti [4J. The CTEQ group has 
studied the magnitude and shape of the strange and antistrange distributions in its global 
analysis |5j, in which the parton distribution function (PDF) set that provided the best 
fit to the global hard scattering data (CTEQ6.5S0) was presented, along with alternative 
sets consistent with the data. The differences for the strange distribution between CTEQ's 
best PDF set and their alternative sets are at the 30% level over the range from x ~ 0.01 
to X ~ 0.1. CTEQ has now published updated PDF sets, CTEQ6.6 [6J, which are an 
improvement over the CTEQ6.5 sets and more suitable for use in the small x region, although 
the physics inputs and assumptions about the strange sea distributions are largely preserved. 
Updated PDF sets have also been published by the MSTW [7] and NNPDF [8J collaborations. 

There are two mechanisms responsible for the quark sea production at experimental 
scales: (I) gluons splitting into quark-antiquark pairs, and (II) non-perturbative contribu- 
tions such as those from the meson-baryon Fock components of the nucleon wavefunction. 
If the masses of m, d and s quarks are light i.e. is much smaller than any relevant ex- 
perimental scale, the sea distributions generated through mechanism (I) can be assumed to 
be flavour independent (SU(3) flavour symmetric) d{x) = u{x) = s{x) and quark-antiquark 
symmetric q{x) = q{x). On the other hand, sea distributions generated via mechanism (II) 
can violate these two symmetries. Mechanism (II) provides a natural explanation for the 
observed SU(2) flavour asymmetry among the sea distributions, i.e. d{x) ^ u{x) [T2|, [13], 
and predicts a strange-antistrange asymmetry p^HTTj . 

We note that, while the strange quark is generally believed to be signiflcantly heavier 
than either u or d {iris ~ 0.1 GeV [H]), this is light compared to experimental scales, and is 
generally ignored in NLO analysis of parton distributions. It is also possible that the larger 
mass of the strange quark could break SU(3) flavour symmetry by restricting the phase 
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space available for ss pairs. This phase space volume is proportional to (Q^ — 4m^), which 
will restrict the phase space when the virtuality of the gluon is close to [12]. However, 
this scale is much smaller than both the typical experimental scales and the initial scales 
used in either approach to fitting PDFs. Even in a model of the nucleon, such as the MIT 
bag model or chiral quark soliton model, where PDFs can be calculated at the model scale 
and evolved up to experimental scales [20], the model scales are typically 0.2 ~ 0.5 GeV^, 
which is still significantly larger than m^. Hence we conclude that mechanism (I), quark pair 
production from gluons, will not give a large contribution to the breaking of SU(3) flavour 
symmetry over the range of scales for which a perturbative approach is reasonable. 

Thus assuming SU(3) flavour symmetry and quark-antiquark symmetry for the sea dis- 
tributions generated via mechanism (I), we can construct the distribution 

A{x) = d{x) + u{x) — s{x) — s{x), (2) 

which has a leading contribution from mechanism (II), and can be calculated using non- 
perturbative models describing that mechanism. A(a;) has been previously calculated in 
the meson cloud model (MCM) [211 122] using covariant perturbation theory and considering 
only baryon-meson components involving pseudoscalar mesons [^ [21]. We present here 
an updated calculation of A(x) using time-ordered perturbation theory in the infinite mo- 
mentum frame, which avoids ambiguities around the use of off-mass-shell structure functions 
[25] . and considering Fock states involving mesons in the pseudoscalar and vector octets plus 
baryons in the octet and decuplet. It is not clear in the MCM where the Fock expansion of 
the nucleon wavefunction should be truncated, as a rapidity gap to identify the final state 
is desirable, however, high mass Fock states are naturally suppressed and only make small 
corrections to the calculation of unpolarized parton distributions. 

Combining our calculation for A{x) with results for the non-strange light antiquark sea 
distributions from global PDF fits we can calculate the total strange distribution S^{x) = 
s{x) + s{x) and the strange sea suppression factor r{x). 

This paper is organized as follows: the theoretical formalism used to evaluate A(x), the 
meson cloud model, is presented in Section 2, the numerical results for A(x), the total 
strange distribution S~^{x), and the strange sea suppression factor r(x) are presented in 
Section 3, and Section 4 is our summary and discussions 
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II. FORMALISM 



Virtual meson-baryon components are created and annihilated continuously in the nu- 
cleon due to energy uncertainty of the system. They are relatively long lived, which enables 
them to contribute to hard processes as these typically have short interaction scales |2T] . 
The presence of the meson-baryon components plays an important role in explaining many 
interesting experimental results, including these from polarized and unpolarized DIS and 
Drell-Yan experiments [I3|. The wave function for the physical nucleon can be written as 

|iV)physical = V^|iV)bare + EE / dy d^k^ ^^^^y, k^) 

BM AA' 

\B\y,k^y,M''{l-y,-k^)) (3) 



where the first term is for a 'bare' nucleon, Z is the wave function renormalization constant, 
and (pBAfiy, k^) is the wave function of the Fock state containing a baryon (B) with longi- 
tudinal momentum fraction y, transverse momentum k^, and helicity A, and a meson (M) 
with momentum fraction 1 — y, transverse momentum — k^, and helicity A'. The probability 
of finding a baryon with momentum fraction y (also known as fluctuation function in the 
literature) can be calculated from the wave function ^'sifiy, k"^ ^ 



fBM/N{y)=Y. dkl<p'A,{y,eMM{y,kl). (4) 

A A' •'^ 

The probability of finding a meson with momentum fraction y is given by 

fMB/N{y) = fBM/Ni^-y)- (5) 

The wave functions and thereby the fluctuation functions can be derived from effective 
meson-nucleon Lagrangians employing time-ordered perturbation theory in the infinite mo- 
mentum frame |25j . 

The mesons and baryons contribute to a hard scattering process such as DIS via the 
Sullivan process |21j, where the virtual photon scatters off a meson or a baryon. The con- 
tribution from the Sullivan processes to the parton distribution functions can be calculated 
via the convolution formula [22] , 

x6q{x) = dyfBM/Nijj) Qb 

= fsM/N^QB (6) 

x5q{x) = /mb/n^Qm (7) 



where qb and qm are valence parton distributions in the baryon and meson. Eqs. ^ and 
([7]) represent the contribution from the process where the baryon and meson respectively 
participate in the hard process (while the other partner in the baryon-meson component is 
a spectator). 

In a recent paper, Strikman and Weiss [26] investigated the dependence on impact param- 
eter of such baryon-meson cloud contributions to sea distributions. They found that about 
one third of the SU(2) symmetry breaking distribution d{x) —u{x) arises from large distance 
contributions, and a similar fraction of A{x) is due to scattering at large impact parameter, 
where they define large distance as greater than the nucleon's transverse axial charge radius, 
estimated at 0.55 fm. We would argue that this approach is contrary to models of nucleon 
structure such as the Cloudy Bag Model (CBM) [27J and the chiral quark soliton model 
xQSM [28j, where quarks couple to pions and other mesons throughout the volume of the 
nucleon, not just near the surface. Also the Sullivan mechanism distinguishes states by their 
rapidity not impact parameter, and it is not surprising that significant contributions to sea 
distributions can arise from e.g. fast moving mesons close to the centre of the nucleon as well 
as those far from the centre. The rapidity gap necessary for identifying MCM contributions 
can be seen in the cross-section for leading neutron production in DIS on protons at HERA 
[29]. 

The leading contributions to the distribution xA{x) come from non-perturbative processes 
(see Eqs. ^ and ([T])). The Fock states we consider include |A^7r) , \Np) , |c<;A^) , lAvr) , |Ap), 
\AK) , |Ai^*) , \T,K), and \T,K*). We obtain the following expression for xA{x) 

xA{x) = (/tt AT/AT + fnA/N + fpN/N + fpA/N + JujN/n) ® K 

- [fAK/N + fAK*/N) ®S^+ (/si^/AT + hK* /n) ® 

+ [fxA/N + fKT,/N + fK*A/N + /i^*S/Af) ® ■ (8) 

In the above expression we have taken the unpolarized parton distributions for the p, and 
u to be the same as that for the tt {Vt^{x)), for which we use the Gliick-Reya-Schienbein 
paramaterization [30] • Also we have taken the unpolarized parton distributions for the 
K and K* to be equal {s^{x)). This distribution, and the baryon valence distributions 
s^(x), s^{x) are calculated using a variant of the MIT bag model developed by the Adelaide 
group [20] and ourselves [3TH3i] . These distributions are calculated at the model scale of 
pI = 0.23 GeV^ and evolved to = 2.5 GeV^ using the programme provided in [35l 136] . 



The expressions for the fluctuation functions involved are presented in Appendix A. The 
calculations details for the fluctuation functions can be found in [ini I2S1 ED E2] ■ We note 
that all the relevant coupling constants are determined by experiment and SU(3) symmetry, 
and the form factors used to make relevant momentum integrals flnite are well constrained 
by flts to DIS data. 

The non-strange light quark sea distributions are well determined by the global PDF flts 
to all available experimental data. Combining the global flt results for d{x) +u{x) with our 
calculation for A(x) we are able to estimate the strange sea distributions via 

X [s{x) + s{x)] = X d{x) + u{x) p.^ — xA{x), (9) 

and the strange sea suppression factor 

r{x) = 1 - ^(^^ _ (10) 

d{x) + u{x) 



Fit 



III. RESULTS 



The numerical results for xA{x) are shown in Fig. 1. The largest contributions to the 
light antiquark sea [d and u) come from the Fock states |A^7r) , \Np) and lAvr), while the 
Fock states involving K and K* mesons and responsible for the strange sea are of roughly 
equal magnitude. We also show the results when Fock states containing K* mesons are 
omitted in order to provide some estimate of the uncertainties in the model calculation. 
We can see that this has a small effect on the calculation of xA{x), and gives us some 
confldence that the uncertainties in the MCM calculations are under control. For comparison 
we show the HERMES measurements for x{s + s) (which are obtained using leading-order 
analysis) combined with the CTEQ group's leading order PDFs (CTEQ6L) for x{d + u) 
to obtain xA(x) as data points. We also show a similar combination of the NuTeV NLO 
parametrisation of x{s + s) [llj with the CTEQ6M PDF set for x{d + m), and results from 
the CTEQ6.5 [5J, CTEQ6.6M 0, and MSTW2008 [7| PDF sets. All results are given 
at = 2.5 GeV^. The shaded area represents the allowed range for xA{x) consistent 
with the uncertainty in the xS~^{x) distribution given by the CTEQ6.5 PDF set [5]. This 
was determined by applying 90% confldence criteria to the di-muon production data sets 
which required that the momentum fraction carried by the strange sea to be in the range of 
0.018 ^ 0.040. 
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FIG. 1: A comparison of the MCM calculations for xA{x) (the thick solid curve - with K* contri- 
bution; the thin solid curve - without K* contribution) with the 'experimental' distributions based 
on HERMES pj (the data points) and NuTeV [llj (the solid curve with open-circle markers) data. 
The resuhs from CTEQ6.6M [6] (the thick dashed curve), MSTW2008 ffj (the dashed curve) and 
CTEQ6.5 [3] (the shaded area with the middle curve giving the central values) are also shown. 



It can be seen that our calculated xA{x) is significantly smaller than the distribution 
calculated using MSTW2008 PDF set, the 'experimental' distribution based on NuTeV 
data, and the central values of the CTEQ6.5 distributions in the region x < 0.2, while 
there is reasonably good agreement with the HERMES results except for the region around 
X ~ 0.10. Our calculation agrees well with the xA{x) obtained using the CTEQ6.6 PDF set 
in the region x < 0.2. We note that the two 'experimental' distributions based on HERMES 
and NuTeV data combined with CTEQ parameterisations of x{d + u) do not agree well. 
It would be very interesting to see if this difference remains when the HERMES data is 
analyzed at NLO in QCD. We note that the MCM calculation of xA(x) is independent of 
any global PDF set for the proton. 

The results for the total strange and antistrange distributions are shown in Fig. 2. For 
these calculations we combined our calculation of xA{x) with the {d + u){x) distribution 
from the CTEQ6.6M set. We find that our calculations are in good agreement with the 
CTEQ6.6M set, and with the HERMES data for the region x < 0.07, but are larger than 
those from the MSTW2008 and CTEQ6.5 sets, and the NLO analysis of NuTeV data [TT]. 
Our calculation for x{s + s) becomes negative for x > 0.25 which is unreasonable. 



8 



^ 0.50 

+ 

0.40 
0.30 

0.20 

0.10 

0.00 
0.1 

-0.10 

X 

FIG. 2: The sum of the strange and antistrange quark distributions, shown as x{s + s), from the 
MCM calculations (the thick solid curve), the HERMES measurements [1] (the data points), and 
the next-to- leading order analysis of NuTeV dimuon data (the thin solid curve). The global 
fit results from CTEQ6.6M [6| (the thick dashed curve), MSTW2008 [T] (the dashed curve) and 
CTEQ6.5 [5J (the shaded area with the middle curve giving the central values) are also shown. 

The reason for this could be that our model calculations overestimate xA{x), or that 
X (d{x) + u{x)^ is underestimated in the CTEQ6.6 set, or both. Our calculations are sensi- 
tive to the calculated PDFs we have used for the baryons (A, A and S), and at large x these 
calculations become unreliable owing to the kinematic limitations of the MIT bag model. 
After convolution with the fluctuation functions, this uncertainty can affect the calculated 
xA{x) and strange distributions in the medium and large x regions, especially as these dis- 
tributions are calculated in terms of differences between the model PDFs. We also note that 
omitting the K* states from the calculation makes almost no difference in this region. 

Our calculation of the strangeness suppression r{x) is shown in Fig. 3. Again, we have 
combined the results from HERMES and NuTeV with CTEQ {d + u) sets at the appropriate 
order to find 'experimental' determinations of this ratio. The MCM calculation gives a 
result that agrees reasonably well with HERMES and the CTEQ6.6M set for x < 0.07, 
and is larger than those obtained in the previous global PDF fits, which suggests that the 
suppression of strange sea relative to the non-strange light antiquark sea is not as large as 
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FIG. 3: The relative strength of strange and antistrange distributions to the non-strange hght 
antiquark distributions, r{x): obtained using the MCM calculations for x{s + s) and CTEQ6.6M 
[B] results for x(d + u) (the thick solid curve), the HERMES [Ij measurements for x{s + s) and 
CTEQ6L [6] set for x{d + u) (the data points), and the next-to- leading order analysis of NuTeV 
dimuon data giving x{s + s) [Hj and CTEQ6M [6] set for x{d + u) (the solid curve). The global 
fitting results from CTEQ6.6M [6] (the thick dashed curve), MSTW2008 f7] (the dashed curve) 
and CTEQ6.5 [5j (the shaded area with the middle curve giving the central values) are also shown. 

has been previously believed. 

IV. SUMMARY AND DISCUSSIONS 

For light quarks, non-perturbative processes are the dominant mechanism for the differ- 
ence A(x) between the non-strange antiquark distributions and the strange sea distributions, 
since the perturbative process of gluons splitting into quark and antiquark pairs is both 
flavour symmetric and quark-antiquark symmetric at leading order. We have calculated this 
difference using the meson cloud model and estimated the total strange plus antistrange 
distribution and the strange distribution suppression factor by combining our calculations 
for A (a;) with the non-strange light antiquark distributions determined from global parton 
distribution functions fits. Our calculations for the strange sea distributions agree with the 
HERMES measurements and that calculated using the CTEQ6.6 PDF set, but are larger 
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than the NuTeV data and those obtained using the CTEQ6.5 and MSTW2008 sets. There 
appears to be a large discrepancy between the HERMES and NuTeV measurements in the 
region x < 0.05, however, distributions in this region are sensitive to the details of the ex- 
perimental analysis. In particular, analysis of the HERMES data at NLO would be required 
in order to check the compatibility of the two measurements of the strange sea distributions. 

This calculation could be extended to helicity dependent sea distributions. We note that 
the hehcity dependent strange distributions xAs{x) + xAs{x) has been measured by the 
HERMES group PQ EH] , and found to be somewhat smaller in magnitude than expected on 
the basis of SU{3)f arguments [32] • Previous MCM calculations [33] have been in accord 
with the HERMES data [3H], however a calculation analogous to the calculation presented 
here could give added confidence in these results. One drawback of this method would 
be that the current data on the helicity dependent non-strange sea distributions has much 
larger uncertainties than the helicity independent data which (along with those of the model 
calculation) may well be larger than any small helicity dependent strange sea distribution. 



Appendix A: Fluctuation functions 



The fluctuation function for — )■ BM can be written as 



47r2y(l - y) (m% - ml^, 



where m^v, ""^s, and uim are the mass for the nucleon, baryon and meson, respectively, and 
m\j^ is the invariant mass squared of the baryon-meson Fock state, 

™|M(y,M = i^^^ + i4^. (A2) 

The vertex includes a form factor G 

kr) = Exp [{ml - m|,,)/(2AL)] , (A3) 

with Acut being a cut-off parameter. In this work we take Aq^^ = 0.8 GeV and A^g^ = 1.0 
GeV for fluctuations involving octet and decuplet baryons respectively. These values of 
the cut-offs give a good description of a wide range of DIS distributions involving MCM 
contributions [121 [161 |3lH33] . 
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For the fluctuations — )■ A^vr, Ai^, and TiK 

V\y,kT) = 3g 



AmBrriNy 



(A4) 



where the couphng constants are taken to be Qn-k/n = 13.07, Qkk/n = —13.98, and Qj^k/n 
2.69 [251 [37]. For the fluctuations ^ Avr, 



(A5) 



where the couphng constant /att/tv = 12.43 GeV~^ [SSI EZ]- For the fluctuations N — )■ A^p, 
A^cj, kK\ and S/sT*, 



1 



, 2[^7 + 2/m^(|/-l)]^|/2^2 



AniBmNy^ \ {y - 1)^ (?/ - 1)^ 

, 2[gy-2fmB{y-l)fy''kl 



~\~krp 



gf(mB - mjv)(|/ - 1)?/ 
"^M(y-l) 

f{y + 1) [fcj. - mB^{y - 1) + y ("-^jv^ly - 1) - "-^M^)] 
muiv - 1) 
^ [fcj. + mBiriN^y - 1)^ - m^f'^y] 
rriMiy-l) 

^ finiNV - rriB) [k^ + niB^ + rriN^y'^ - {niB^ + + ^1^"^)] 

muy 



+2y^ 



g{y - l)(mB - mjyy) 

^ 2/ [kl - (mj? + mAr)(?/ - l)(mB - mAft/)] 



n 2' 



(A6) 



(y-l) 

where the couphng constants are taken to be Qnp/n = 3.25, fNp/N / Qnp/n = 6.1/4mjv, 
S'A^w/Af = 10.09, Jnuj/n / Qnui/n = 0, gAK*/N = —5.63, f\K*/N = —4.89 GeV"\ gT.K*/N = 
-3.25, and /si^*/^ = 2-09 GeV-^ j25l|37]. For the fluctuations A^ -> Ap, 

f 



V\yM) 



QmB^rriNiy — 1)^?/ 



-— ^ {fc^ + [(4y2 - 4y + 3) niB^ + (2mM' + m^v') y^ 



+ 



- 1)2 (3?/2 - 2?/ + 3) + 2|/2 [(^2 ^ 2) ^m' + 2mN\y - 1 



+mM rriN y + mB (y - 1) + tub tum [3mM + 2mN {y -I] 
-l2mB^mM^rnN{y - 1) V + m^' [2mM' + 3miv'(l/ - 1)'] (y - ify^] ,(A7) 



12 



where the couphng constant /ap/n = 20.82 GeV ^ j25| |37]. 
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